Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN KTeX style file v2.2) 



On the X-ray spectra of luminous, inhomogeneous accretion flows 



A. Merloni^*, J. Malzac^, A.C. Fabian^ and R.R. Ross^ 

^ Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild-Strasse 1, D-85741, Garching, Germany 
^Centre d'Estude Spatiale des Rayonnements, CNRS-UPS 9, Avenue du Colonel Roche, 31028 Toulouse, Cedex 4, France 
- - - ^Institute of Astronomy, Madingley Road, Canbridge, CB3 OHA 

, "^Physics Department, College of the Holy Cross, Worcester, MA 01610, USA 

o . 

o ■ 

^ . 5 February 2008 

ABSTRACT 

' We discuss the expected X-ray spectral and variability properties of black hole accretion discs 

at high luminosity, under the hypothesis that radiation pressure dominated discs are subject 
to violent clumping instabilities and, as a result, have a highly inhomogeneous two-phase 
structure. After deriving the full accretion disc solutions explicitly in terms of the parameters 
of the model, we study their radiative properties both with a simple two-zones model, treat- 
able analytically, and with radiative transfer simulations which account simultaneously for 
\ energy balance and Comptonisation in the hot phase, together with reflection, reprocessing, 

\^ ■ ionization and thermal balance in the cold phase. We show that, if not only the density, but 

' also the heating rate within these flows is inhomogeneous, then complex reflection-dominated 

i spectra can be obtained for a high enough covering fraction of the cold phase. In general, 

■ large reflection components in the observed X-ray spectra should be associated with strong 

soft excesses, resulting from the combined emission of ionized atomic emission lines. The 
Q^. variability properties of such systems are such that, even when contributing to a large fraction 

I ' of the hard X-ray spectrum, the reflection component is less variable than the power-law like 

O \ emission originating from the hot Comptonising phase, in agreement with what is observed in 

many Narrow Line Seyfert 1 galaxies and bright Seyfert 1 . Our model falls within the family 
^ ' of those trying to explain the complex X-ray spectra of bright AGN with ionized reflection, 

but presents an alternative, specific, physically motivated, geometrical setup for the complex 
multi-phase structure of the inner regions of near-Eddington accretion flows. 
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1 INTRODUCTION 

Black Hole systems, either in X-ray binaries (BHXRB) or in active 
galactic nuclei (AGN), when shining at luminosities close to the 
Eddington limit are thought to be powered by accretion through 
geometrically thin, optically thick discs, where angular momentum 
is transferred by the stresses due to magneto-rotational instability 
(MRI; Balbus & Hawley 1991, 1998). From the observational point 
of view, the presence of such discs is inferred from the spectral 
properties of transient BHXRB at (or close to) the peak of their out- 
burst luminosities (in the so-called High/Very High/Intermediate 
states, see e.g. Belloni 2004, or McClintock & Remillard 2005 for a 
different nomenclature, and references therein) and from the dom- 
inance of the quasi-thermal UV emission in bright Quasars (see 
e.g. Malkan 1983; Laor 1990; Shang et al. 2005). According to the 
standard models (Shakura & Sunyaev, 1973), such accretion discs 
are radiatively efficient, as they are able to convert into radiation a 
sizeable fraction (between 6 and ~40%, depending on black hole 
spin) of the available gravitational energy of the accreting gas. In- 
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deed, radiative efficiency estimates obtained by comparing the total 
supermassive black hole mass density in the local universe with the 
accretion power released by Quasars and AGN suggest that this is 
probably the main mode of growth by accretion of supermassive 
black holes (see e.g. Fabian & Iwasawa 1999; Elvis, Risaliti and 
Zamorani 2002; Yu and Tremaine 2002; Marconi et al. 2004; Mer- 
loni, Rudnick & Di Matteo 2004). It is therefore clear that a good 
understanding of the accretion mechanism associated with nearly- 
Eddington sources can shed light on fundamental issues of black 
holes astrophysics. 

Nonetheless, according to the above mentioned standard ac- 
cretion disc solutions, such highly luminous discs should be radi- 
ation pressure dominated and therefore unstable to perturbations 
of both mass flow (Lightman & Eardley 1974) and heating rates 
(Shakura & Sunyaev 1976). Thus, it is not clear yet to what ex- 
tent this standard solutions represent a realistic description of the 
observed systems. In recent years, both analytic works (Blaes and 
Socrates 2001, 2003) and simulations (Turner, Stone and Sano 
2002, Turner et al. 2003b, Turner 2004) have shown that mag- 
netized, radiation pressure dominated accretion discs may be in 
fact subject to violent clumping instabilities if magnetic pressure 
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exceeds gas pressure and photons are able to diffuse from com- 
pressed regions. Large density variations may also be caused by 
photon bubble instabilities (Gammie 1998), which may develop 
into a series of shocks propagating through the plasma (Begelman 
2001,2002; Turner at al. 2005). These instabilities may in turn have 
profound effects not only on the nature of the cooling mechanism 
of luminous discs, but also on their observational appearance (see 
e.g., Davis et al. 2004; Ballantyne et al. 2004), as we discuss in this 
work. 

Radiative models for inhomogeneous, clumpy discs that ex- 
plain the X-ray properties of accreting black holes have been pro- 
posed since long time. Most of the previous works were devoted to 
the study of the internal state of the cold clumps and on their radia- 
tive output (Guilbert & Rees 1988; Celotti, Fabian & Rees 1992; 
Collin-Souffrin et al. 1996; Kuncic, Celotti & Rees 1997). Krolik 
(1998) studied analytically the overall equilibrium of a clumpy two- 
phase accretion disc, while Malzac & Celotti (2002; MC02) stud- 
ied numerically the Comptonization equilibrium between the cold 
clumps and the hot surrounding plasma and calculated the emitted 
X-ray spectrum. 

In this paper, we explore in detail the consequences of the hy- 
pothesis that accretion flows close to the Eddington rate are indeed 
inhomogeneous. In particular, we present a detailed study of their 
expected X-ray spectral and variability properties, extending the re- 
sults of MC02 both numerically and analytically. From the numer- 
ical point of view, we calculate the emerging spectrum from an in- 
homogeneous, two-phase, clumpy flow by coupling the non-linear 
Monte Carlo code of Malzac & Jourdain (2000), that calculates the 
Compton equilibrium between the cold and the hot phase, with the 
X-ray ionized reflection code of Ross & Fabian (1993, 2005), that 
accurately computes the reprocessed radiation in the cold phase. 
From the analytic point of view, we are able to show that a crucial 
factor in determining the broad band spectral properties of an in- 
homogeneous flow is not only the amount of cold clouds/filaments 
pervading the hot plasma (as considered in all the previous studies 
cited above), but also the level of inhomogeneity of the heating in 
the hot phase itself. Indeed, we will demonstrate both analytically 
and numerically, the very important result that the more localized 
the heating is (as due to, e.g. a reconnection event), the more the 
emerging spectra are dominated by reflection/reprocessing, as al- 
ready suggested by Fabian et al. (2002). 

In this paper we will also show that a model for an accretion 
flow which is inhomogeneous both in its density and heating struc- 
ture can explain many of the observed properties of black holes 
accreting at rates close to the Eddington one (e.g. bright QSOs 
and Narrow Line Seyfert galaxies). In particular, flows with vary- 
ing degrees of inhomogeneity naturally produce little variability in 
the observed reflected component associated with large variations 
of the power-law continuum, and reflection dominated spectra at 
the lowest-luminosity levels, in a strikingly similar fashion to what 
predicted by the so-called light-bending model (Miniutti & Fabian 
2004). Similarly, such a configuration naturally produces multi- 
ple reflections of the X-ray emission that can enhance the broad 
emission and absorption features in the observed spectra (Collin- 
Souffrin et al. 1996; Ross, Fabian & Ballantyne, 2002). 

The paper is structured as follows: In Section|2|we first intro- 
duce our simple toy model for a clumpy, two-phase medium with 
inhomogeneous heating. As this model is treatable analytically, we 
are able to discuss the main expected properties (spectra and vari- 
ability) over a wide area of the parameter space. In Section|3|we 
then discuss our numerical simulations, while in Section|4|we dis- 
cuss the relevance of our results for observations of high accretion 




Figure 1. A simple two-zone model for an inhomogeneous accretion flow. 
The hot plasma emitting the hai'd X-rays is pervaded by small clouds of 
dense matter, optically thick to Thomson scattering. Both cold clumps and 
the hot plasma are uniformly distributed, but the heating is only concen- 
trated inside the inner sphere of radius Z = hH. The system is assumed to 
be in radiative equilibrium. 

rate BHXRB and AGN and their implications for models of ac- 
cretion flows. Finally, we draw our conclusions in Section |5| In 
addition, we have developed in Appendix |3 a full inhomogeneous 
accretion disc model based on the picture discussed in Section |2| 
in which the dependence of our model parameters on radius, black 
hole mass and accretion rate is made explicit. 



2 A SIMPLE TOY MODEL 

In order to gain insight into the main expected observable proper- 
ties of inhomogeneous discs, we introduce first a simple, yet moti- 
vated, idealized model, which is tractable analytically. In Appen- 
dix 1^ instead, we present a full inhomogeneous accretion disc 
model, closely following the ideas of Krolik (1998), with the aim 
of constraining the values and the scalings of the parameters we 
introduce in order to describe the disc inhomogeneity (see below). 

Let us consider a sphere of hot plasma with uniform density 
and radius H, pervaded by small lumps of cold matter randomly 
distributed throughout the hot phase. Such a system may represent 
the inner region of an accretion flow (with H ~ few Schwarzschild 
radii, iig), or a patch of the accretion disc itself in the radiation 
pressure dominated region, where clumping instability is at work 
(with H ~ disc scaleheight, see Appendix |A}. Then we can write 
H — rRs = 3 X W^^mrro cm, where mr is the mass of the central 
black hole in units of IO^Mq, and ro = r/l} As already discussed 
in MC02, the lack of strong absorption features in the X-ray spectra 
of accreting black holes suggests that any putative cold clump of 
matter immersed in the hot plasma should be optically (Thomson) 
thick (see also, Kuncic et al. 1997), with zero transmission. Thus 
we consider clouds with column densities > 10^^ cm~^, which 
implies a density rici > 3.3 x 10^''(m7e_iro)~^ cm""^, where we 
have introduced the parameter e — 0.1e_i, the ratio of the cloud 
size to the system size H. 

We assume that the whole system is in radiative equilibrium. 

^ Throughout the paper we use the representation Ax = A/W^. 
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Figure 2. The relative intensity of the emerging hard X-ray Comptonized 
(LCi upper left panel), reprocessed UV (Luv, lower left panel), and re- 
flected (Lp , upper right panel) luminosities is plotted as color-coded con- 
tours in the (tb , h) plane. Superimposed on them, in the two upper panels 
are the contours of constant reflection fraction (solid hnes), ranging from 
= 0.5 to = 4, while in the lower panel we plot the contours (dashed 
lines) of constant power-law indexes F, from 1.9 to 2.4. The optical depth 
of the hot phase is assumed to be unity r-p = 1 



However, differently from MC02, we allow the heating to be lo- 
calized in a small region, rather than diffused throughout the hot 
phase. For simplicity, we assume that the heating of the hot plasma 
takes place only in an inner concentric sphere of radius Z = hH 
(see Fig. 0, where a power ih is deposited. The outer zone can 
therefore be viewed as a passive reprocessor (see Appendix 151. In 
our idealized model the injected power eventually escapes the sys- 
tem as three different kinds of radiation: 



Lh = I/T 



Lc 



(1) 



respectively the UV/soft X-ray reprocessed luminosity, the Comp- 
tonized (power-law like) luminosity and the reflected luminosity. 
In our calculation we also assume that the UV/Soft X-ray radia- 
tion reprocessed by the clouds is the only source of soft photons 
(no internal dissipation in the cold clumps, see e.g the discussion in 
Collin-Souffrin et al. 1996, or Malzac 2001). 

If the clouds are sufficiently small, we can use an effective 
cloud optical depth, tb, as the main parameter controlling the 
amount of cold matter. This is defined in such a way that a photon 
crossing the hot plasma has a probability 1 — e^^^ of intercept- 
ing a cloud. The second parameter that affects the emerging spec- 
trum is the optical depth to Thomson scattering of the hot plasma, 
tt. Comptonization models for the X-ray emission from compact 
objects typically constrain the hot plasma Thomson optical depth 
Tt ~ 1 (see e.g. Zdziarski 1999). This is also consistent with the 
upper and lower limits on tt derived by Krolik (1998) from the 
thermal and evaporative equilibrium of a clumpy disc. 

A rough order of magnitude estimate of the expected ioniza- 
tion parameter can be derived as follows: radiatively heated clouds 
with no internal dissipation, if dense enough (as we are assum- 
ing here, see also Collin-Souffrin et al. 1996; Kuncic et al. 1997), 



will re-emit the absorbed flux as a black-body with typical energy 
fcTbb — 30Z/4{^(m7ro)~^'''^ eV, where L is the bolometric lumi- 
nosity of the system and the weak dependence on the albedo has 
been neglected. Then, assuming that in the hot phase gas pressure 
dominates over radiation pressure (an assumption that can be veri- 
fied a posteriori, see details in Appendix]^, thermal pressure bal- 
ance between the cold clouds and the hot phase implies that 
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Therefore, the typical ionization parameter^ 

C = 4.i./n. . 3.5xlO^L^/^^..o)-^^ {^f 

It is worth stressing, however, that in the inhomogeneous heating 
model we are discussing, a range of ionization parameters is ex- 
pected, given the sensitivity of ^ to the actual geometry of the sys- 
tem. The above equation has to be considered just an order of mag- 
nitude estimate. Indeed, the above estimates are consistent with the 
values obtained with the full disc modeling outlined in AppendixlAl 
where the explicit scalings of cloud density, temperature and ion- 
ization parameters with accretion disc radius, black mass and ac- 
cretion rate are presented. In the following (see Section|3}, we will 
explore the effects of changing the value of ^ on the spectra emerg- 
ing from our inhomogeneous flow. 

The other parameters needed to compute the emerging radia- 
tion intensity are the energy and angle-integrated single-reflection 
albedos for a Comptonized and reflection spectrum, a and or, re- 
spectively. Malzac, Dumont & Mouchet (2005) have shown that, 
when the reprocessed spectrum is computed self-consistently the 
X-ray albedo is a strong function of the ionization parameter and 
of the spectrum of the illuminating radiation. Moreover, multiple 
reflections start to become important for tb ~ 1. In a simple one- 
zone model, the net albedo produced by multiple reflections should 
be approximately given by (Eq. 2 of Ross, Fabian and Ballantyne 
2002): 



anct = ea/[l — (1 — e)a], 



(4) 



where e is the escape probability for a reflected photon and a is 
the albedo for a single reflection by the material. We should expect 
that, whenever heating is inhomogeneous the cold clumps in the 
inner heated zone and in the outer one will not only be exposed 
to illuminating radiation with different spectra, but also experience 
different numbers of multiple reflections. Thus, it is plausible to 
conceive that in general e will not be identical to exp(— tb) as 
would be in a simple one zone model. For the sake of simplicity, 
and because the main purpose of the toy model discussed here is 
to give a coherent qualitative picture of the properties of inhomo- 
geneous flows, we have chosen to calibrate the net albedo directly 
on the simulated spectra that will be presented in Section |3| To a 
very high degree of accuracy, we find that, for a ionization param- 
eter of ^ = 300, the net albedo is a function of the optical depth of 
the cold clumps only, and can be expressed with a similar expres- 
sion by replacing e = exp(— rB/1.43), with a — 0.36 (see also 
figure^. The reflected radiation can itself be reflected many times 
for the most extremely inhomogeneous discs. Also for its albedo 
we have used the expression J4j with e = exp(— tb/1.82), with 
Or = 0.4. This latter value is also consistent with the Monte Carlo 



^ Here, as in the remaining of the paper, the ionization parameter will be 
expressed in erg cm s~ ^ . 
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simulations of Malzac (2001), where the albedo for a typical reflec- 
tion spectrum illuminating neutral matter was computed. We have 
also verified that assuming an increased albedo of or = 0.6 to 
account for the effects of ionization have negligible effects on the 
final results of our analytic two-zones model. 

We can then proceed as in Section 4 of MC02 and solve for 
the three different luminosities Luv, Lc and Lr (and the corre- 
sponding reflection fraction R) making use of the analytic approxi- 
mations for the corresponding photon escape probability calculated 
from Monte Carlo simulations (see Eqs. A3 and A5 in MC02). We 
keep the heating rate, Lh, fixed and explore the parameter space 
spanned by 0.1 < tr < 8 and 0.05 < h < 0.9. Details of the 
analytic calculation of the emerging luminosities in this two-zone 
model are given in the Appendix IHl Once the soft and the Comp- 
tonized luminosities are found, we determine the spectral shape 
of the X-ray emission by calibrating the power-law index F with 
the numerical simulations (see below), where the power law index 
of the Comptonized emission is calculated by fitting the emergent 
spectra with the PEXRAV model in the 2-30 keV range. Apart from 
a very weak dependence on the heating inhomogeneity parameter, 
h, the power-law index depends only logarithmically on re, with 
the spectrum being softer for larger cold clouds column densities. 
From the analysis of the results of the simulations, we derive the 
following phenomenological expression for the power-law index 
as a function of tr : F ~ 2 + 0.4 log tr . 

Figure |2| shows the emerging luminosities as coloured con- 
tours in the two-dimensional parameter space (tr, h), for a Thom- 
son optical depth of the hot phase tt = 1. Shown are the relative 
intensities (normalized to the total heating rate Lh) of the X-ray 
Comptonized luminosity (Lc, upper left panels), of the reflected 
luminosity (Lr, upper right panels) and of the reprocessed UV/soft 
X-ray one (Luv, lower left panels). Superimposed on them, in the 
two upper panels are the contours of constant reflection fraction, 
defined R = Ln/{anctLx), where Lx = Lc{E > IkeV) is 
the Comptonized luminosity emerging above 1 keV. This has been 
calculated analytically from Lc assuming that the Comptonized 
emission is power-law with index F extending from fcTbb = 50 eV 
up to an energy corresponding to the temperature of the hot elec- 
trons e = fcThot/(moC^) ~ (3/19)/[(F - 1)t^^^] (Wardzinski 
and Zdziarski 2000). In the lower panel we plot the contours of 
constant power-law spectral index, F. 

As expected, the reflection fraction 7? increases with tr, as 
does the relative intensity of the soft reprocessed emission: the 
larger the filling factor of the cold component, the more reflection 
dominated the emerging spectrum will be. In general, large reflec- 
tion fractions are always expected when the reprocessed component 
is strong compared to the direct Comptonized X-ray emission. 

The relative intensity of the soft reprocessed luminosity, Luv, 
is larger in the centrally heated case than in the uniformly heated 
case, if tr is large. As a result, for large enough tr and localized 
heating (i.e. in the upper left comer of the plots shown in Fig.|2j, the 
X-ray continuum power-law is dominated by the reflection compo- 
nent and by the thermal reprocessed emission. 



2.1 On the variability of the reflection component 

One more interesting fact emerges from the study of this simple 
two-zone toy model. As can be seen in Fig. |2| the reflected lu- 
minosity Lr does not vary as much as Lx and Luv across the 
explored region of the parameter space, at least as long as tr is 
smaller than ~ 2 (note that the color scale of Fig.|2|is linear). In 
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Figure 3. The relative intensity of the reflection component (/?, solid line), 
of the reprocessed UV (Luv, dot-dashed line) and reflected (Lr , dashed 
line) luminosities are plotted as functions of the emergent X-ray (Comp- 
tonised power-law) luminosity above 1 keV (Lx) for a varying cloud op- 
tical depth Tr . All luminosities are renormalized by the total heating rate 
L^, which is here assumed to be concentrated in an inner sphere of radius 
/i = 0.1. The vertical dotted lines mark the X-ray luminosities correspond- 
ing to TB =0.5, 1,3,5. 



fact, the reflected luminosity Lr has a broad maximum almost co- 
incident with the locus of J? ~ 1. To make this clearer, we have 
plotted in Fig.|3|the reflected luminosity, the soft reprocessed lumi- 
nosity and the reflection fraction as a function of the X-ray (Comp- 
tonized) luminosity above 1 keV for the case of inhomogeneous 
heating /i = 0.1 and tt = 1 (those are also the values of the pa- 
rameters adopted in our numerical simulations of Section|3|below). 
The fact that the reflected luminosity has a maximum, implies that 
large variations of the emergent X-ray luminosity, Lx, associated 
with changes in the cold clump integrated optical depth correspond 
to only modest variations of the reflection component, at least as 
long as Lx/Lh > 0.05. On the other hand, for low values of the 
Comptonized X-ray luminosity, the reflected luminosity correlates 
with Lx, while at high Lx, the trend is the opposite. 

This global behaviour of Lr vs. Lx is analogous to that of the 
so-called "light bending" model (Miniutti & Fabian 2004), which 
has been previously invoked to explain the variability properties of 
the continuum and of the iron line in a number of Narrow Line 
Seyfert 1 (NLSl) galaxies and galactic black holes in the very high 
state (see, e.g. Miniutti et al. 2003; Miniutti, Fabian & Miller 2004; 
Fabian et al. 2004; Rossi et al. 2005). In particular, the defining 
property of the reflected luminosity to show a broad maximum as 
a function of the emerging X-ray (Comptonized) one is a common 
characteristic of both models. We will return to this similarity later 
on. 

We have explored the variability property of our simple toy 
model in greater detail. In order to do so, we have computed 
the fractional variability cr; of each spectral component (i — 
X, UV, R), associated with variations 5 log tr , 5 log h of the cold 
clouds optical depth and of the inner heated region size, respec- 
tively: 



0-45 log TR , 5 log /l) = 



d log Li 



5 log TB + 



d log Li 



I (5 log ft. (5) 



d log TR d log h 

Figure|4|shows as contour plots the fractional variability levels 
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Figure 4. The total variability fraction, associated with variations of 10% 
of TB and h = Z/H, of Comptonized X-rays emission above 1 keV (crx, 
upper left panel), reprocessed UV (o"uv> lower left panel), and reflected 
(iTj^, upper right panel) luminosities is plotted as color-coded contours in 
the (tb, h) plane. Superimposed on them, in the two upper panels are the 
contours of constant reflection fraction (solid Unes), ranging from R = 0.5 
to = 4, while in the lower panel we plot the contours (dashed lines) of 
constant power-law indexes F, from 1.9 to 2.4, which are the same as in 
Fig.H 

of the three components calculated for 5 log t-b = 5 log h — 10%. 
The fractional variability level of the X-ray emission can thus be 
relatively high (up to ~ 50% for tb ~ 2) even for relatively 
small variations of the total optical depth. According to our sim- 
ple toy model of spherical clouds in a spherical volume, the total 
number of clouds, A^, is related to the cold phase optical depth by 
Tb = (3/4) A'^e^ . The total number of clouds will have a dispersion 
around its mean value of ~ \/7V, so that Atb/tb — 0.1 corre- 
sponds to A'' ~ 100 and e ~ O-II-^tb- Interestingly, one obtains a 
similar constrain on the cloud size from considering thermal equi- 
librium in the presence of heat conduction (see Krolik 1998, eq. 5): 
e ~ 0.2r~'(fcThot/100 keV)^''^ 

Overall, the fractional variability of Lx is higher if the reflec- 
tion fraction R is high. For tb < 2, and R < 2, the variability of the 
Comptonised emission, crx, is more than twice that of the UV/soft 
X-ray and reflected components. For even lower values of the cold 
phase optical depth tb < 0.5, and correspondingly lower reflection 
fractions, the amount of variability is comparable in both hard and 
soft, reprocessed, components, while for the opposite extreme case 
of very large tb (i.e. for the case in which most of the X-ray and 
reflected emission is hidden from the view), the reflection compo- 
nent at high energies will vary most dramatically, while the UV/soft 
X-ray emission should be almost constant, with the Comptonized 
power-law emission varying by as much as 60-70%'^ 

The direct analogy with the light bending model presented 

Obviously, it should always be kept in mind that we are here considering 
only the variability induced by structural/geometrical changes at fixed total 
heating rate Ljj . One should naturally expect, in a realistic situation, that 
the heating rate is also variable on a dynamical timescale, and its variations 
should contribute to the observed overall variability. 



above lends itself to a simple geometrical explanation. In the more 
general framework of two-phase models for the X-ray spectra of 
accreting black holes, the main spectral and variability properties 
are determined by the radiative feedback between the hot and the 
cold phases. Such feedback is strongly dependent on the geometry 
(and on the topology) of the two phases and, in particular, on the 
sky covering fraction of the cold phase as seen by the hot, Comp- 
tonising medium. Reflection dominated spectra are expected when 
the cold phase intercepts most of the photons coming from the hot 
phase. This, in the light bending model, is achieved via general rel- 
ativistic (GR) effects, while here it is a result of the clumpy and 
inhomogeneous nature of the inner disc. 

In more general terms, one can define as "open" any geome- 
try in which the hot. X-ray emitting plasma is photon starved (i.e. 
patchy coronae, inner ADAF -I- outer discs etc.). Such a geometry 
will produce hard X-ray spectra, little soft thermal emission and 
weak reflection component (see e.g. the simple model of Zdziarski, 
Lubiiiski & Smith 1999). On the other hand, a "closed" geome- 
try will instead correspond to a very large covering fraction of the 
cold phase, with associated strong soft emission, softer spectra and 
strong reflection fraction (Collin-Souffrin et al. 1996). Such a ge- 
ometry corresponds, observationally, to the highly accreting black 
holes. Whether this is due predominantly to strong GR effects as 
expected from centrally concentrated coronae atop a cold, thin disc 
in the immediate vicinity of the black event horizon, or to the more 
mundane inhomogeneous structure of a radiation pressure domi- 
nated disc, remains to be tested with accurate spectral modeling. 

Is it possible to discriminate between these two scenarios for a 
"closed" geometry? In principle, the relativistic blurring induced by 
the differential rotation of the inner disc should always be taken into 
account when fitting observed spectra. In the original light-bending 
model, where the illuminating source is a point-like source above a 
standard geometrically thin disc, the ratio of the reflected compo- 
nent to the power-law continuum is determined by the same effects 
that determine the shape of the relativistic lines, while if the disc is 
truly inhomogeneous, the two effects can be decoupled. Therefore, 
simultaneous spectroscopic studies of relativistically blurred emis- 
sion lines and of the broad band continuum and variability could 
be effectively used to discriminated between a pure light bending 
model and a clumpy disc. Detailed predictions for the latter, how- 
ever, require the combination of sophisticated MHD and radiative 
transfer simulations. 



3 NUMERICAL SIMULATIONS 

The toy model presented in Sec. |2| provides important qualitative 
predictions regarding the strength and variability of the different 
spectral components emerging from a clumpy accretion flow. In 
this Section, we perform detailed radiative transfer calculations to 
investigate quantitatively the shape of the predicted spectra. These 
calculations account simultaneously for energy balance and Comp- 
tonisation in the hot phase, together with reflection, reprocessing, 
ionization and thermal balance in the cold phase. The radiative pro- 
cesses in the hot gas are simulated using the Non-Linear Monte 
Carlo Code of Malzac & Jourdain (2000), the reprocessed radia- 
tion from the cold clouds is computed using the photo-ionisation 
code of Ross & Fabian (1993,2005). The results of each code are 
used as input for the other code and self-consistency is achieved 
by multiple iterations as described in Malzac, Dumont & Mouchet 
(2005). 
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H 

Figure S. In numerical simulations, the disc is modeled as an infinite slab 
made of a paving of identical cubic regions. As shown in this sketch, each 
cube is itself divided into 10 cells to account for inhomogeneity. The cube 
is uniformly filled with hot/ionized plasma and cold clumps, but dissipation 
occurs only in the central (colored) zone. 



3.1 Setup 

The geometry we consider for the purpose of our numerical simu- 
lations differs from that of the simple analytical model presented in 
Section|2| The basic idea behind our calculations is to study the ob- 
served spectra emerging from a clumpy accretion flow in which the 
heating is localized. Thus, the schematic two-zone toy model dis- 
cussed in Section 2 represents here a patch of the accretion disc, 
with H ~ disc scaleheight (spherical symmetry and concentric 
zones were assumed there only to allow an exact analytic treat- 
ment of the two zone model). In the simulations, the accretion disc 
is modeled as an infinite slab made of a paving of identical cubes of 
size H and dissipation takes place only close the center the cubes. 
The hot plasma density is uniform and characterized by a Thomson 
optical depth tt = UeCttH. The disc is also uniformly pervaded 
by cold dense clumps with an effective cloud optical depth tb. In 
the case of spherical clumps the total number of clumps in the cube 
A'^ is given by re = ivNe^ (see MC02). In this idealization, we 
can afford to simulate the emission from only one paving stone by 
assuming that the same radiation flux that escapes to a neighboring 
cube through one side, re-enters on the opposite side. 

To account for expected gradients of temperature of the hot 
plasma and ionization parameter of the cold clouds, the cube is di- 
vided into ten cells. Those cells are defined as the intersection be- 
tween 5 equidistant planes that are parallel to the accretion disc and 
a cylinder the axis of which is normal to the slab-disc and crosses 
the center of the cube (see Fig.|5j- The plasma is heated only in 
the central cylindrical cell. The temperature in each cell is com- 
puted from the local energy balance. For a given volume averaged, 
ionisation parameter 5, the ionisation gradient is computed from 
the local radiative energy density. The compactness of the active 
region is parametrized through h = Z/H, where Z is the radius 



of the cylinder. We assume that the cloud have a constant density, 
rici — 10^^ cm"'', and standard (solar) abundances (see Ross and 
Fabian 2005, and references therein for the atomic data used in the 
code). The values for the cold clouds density and ionization pa- 
rameters are therefore chosen within the range of values expected 
for accretion flows around supermassive black holes (see Appen- 
dix A). The explicit scaling with mass, and the possible relevance 
of our model for stellar mass black holes are discussed in 14.31 

3.2 Results 

Fig.|6|shows a sample of model spectra that illustrates the effects of 
varying tb from 1 to 8. The emerging spectra are clearly very com- 
plex, with prominent signatures of ionized reflection dominated by 
a complex of UV and soft X-ray emission lines. Strong ionization 
edges are observed for low ionization parameter (for high values of 
^ Comptonization has a relatively large influence on the emerging 
spectra and help smearing out ionization edges and lines). 

As mentioned earlier, the complex geometry of our setup will 
inevitably lead to multiple reflections, which in turn reduce the net 
albedo of the cold clumps. We have studied in detail this effect by 
calculating the total albedo as the ratio of incident (Comptonized) 
flux to the reprocessed one at energies above 0.1 and 1 keV. The 
results are shown in Figure^ 

In order to quantify at least the most prominent characteristics 
of these spectra in the X-ray band, we have fitted them in the 2-30 
keV range under XSPEC with the PEXRAV model plus a Gaussian 
line. During the fit the high energy cut-off energy Ec and the in- 
clination angle were fixed at 400 keV and 30 degrees respectively, 
and the abundances were fixed to solar. Obviously, we are not in- 
terested here in the goodness of these fits, rather we want to use the 
fit parameters as measures of the spectral complexity. The results 
of the fitting procedure confirms that reflection dominated spec- 
tra (in the sense usually implied by an observer, i.e. spectra with 
R > 1, where R is the value of the reflection fraction parameter 
of the PEXRAV model) can be produced with our model for the 
inhomogeneous disc structure. For tb = 8 the reflection fraction 
is R ~4-5, even at large ionisation parameters (^ ~ 3000). When 
comparing these results with the simple toy model of Section|2| it is 
apparent that the PEXRAV fits indicate somewhat larger reflection 
fractions at small tb ■ 

Below ^ 1 keV, all spectra show a distinct curvature that, 
when plotted as a ratio to the best fit power-law in the 2-10 keV 
range, would appear as a so-called "soft excess", indeed a ubiqui- 
tous feature in the spectra of QSOs and Seyfert galaxies (see e.g. 
Wilkes & Elvis 1987; Piconcelli et al. 2005). Within our model, 
such emission is accounted for by a plethora of emission lines pro- 
duced by ionized reflection of the Comptonised continuum off the 
cold clumps as proposed earlier by Ballantyne, Iwasawa and Fabian 
(2001), Fabian et al. (2002), and more recently by Crummy et al. 
(2005). 

Figure |8| shows in red a detailed view of the emerging spec- 
trum for the case tb = 3, ^ = 300. By looking at the ionization 
structure within the gas, it is apparent that the gas is still highly 
ionized at the outer surface, with the dominant ions being C VII, 
N VIII, O IX, Ne X, Mg XI, Si XII, S XIV and Fe XVIII. How- 
ever, the ionization decreases rapidly with depth, and by a Thomson 
depth of 1 .4, the least-ionized species treated dominates for each el- 
ement. Specifically they are, with their ionization potentials: C III 
(47.9 eV); N III (47.5 eV); O III (54.9 eV); Ne III (63.5 eV); Mg 
III (80.1 eV); Si IV (45.1 eV); S IV (47.3 eV) and Fe VI (99.1 eV). 
Within our code, the lowest absorption edge treated is thus at 45 
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Figure 6. Effects of tb and ^ on the emerging spectrum. Angle averaged spectra calculated for h = 0.1 and tt = 1. In each panel, are shown the 
corresponding tb and ^ as well as the best fit parameter R and F obtained when these spectra are fitted with pexrav in the 2-30 keV range and the optical 
to X-ray spectral slope Oox- The thin dashed and soUd curves indicate the Comptonised and the reprocessed spectra respectively. The total spectra are shown 
by the thick solid curve. 
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f=300, T^=l 




Figure 7. Solid symbols show the net albedo of the cold clumps calcu- 
lated from our simulations, after multiple reflections are taken into account. 
With triangles we show the results obtained by considering only the repro- 
cessed radiation above 1 keV, while squares show the results obtained by 
considering reprocessed radiation above 0.1 keV. Solid and dashed lines 
show the approximations to the analytical results by the analytic expres- 
sion J4j, with a = 0.36 (equal in both cases), e = exp(— tb /1.43) and 
e = exp(— tb/1.82), respectively. The solid line also shows the relation 
used in the analytic two-zone toy model described in Sectionl2l 



eV. Below that energy, the only absorption treated is free-free ab- 
sorption, and radiation can enter and leave the gas with ease. We 
should note here that in the code helium is assumed to be fully ion- 
ized (so that its only contribution is to the free electron density). 
However, since the least-ionized species of C, O and Si appear to 
be present in this particular model, some He II might be present in 
reality, which would produce an edge at 54.4 eV if it were included 
in the calculation. The four strongest low-energy emission lines are 
due to Si IV (1 1.0 eV), O III (17.6 eV), C III (21.6 eV), and O III 
(33.1 eV). 

The spectra shown in Figure|6|are those emerging from a sys- 
tem at rest. It is obvious that, when considering the application of 
such a model to the inner region of an accretion disc, the effects 
of both the turbulent motion of the clumps and the overall rotation 
pattern of the two-phase disc should be taken into account (as well 
as cosmological redshift, obviously). Indeed, in all previous direct 
applications of ionized reflection models to observed data (see e.g. 
Fabian et al. 2002; Crummy et al. 2005), the absence of sharp emis- 
sion lines both in the soft X-ray and in the Iron Kq line regions of 
the observed spectra was taken as an indication of extreme rela- 
tivistic blurring in the inner region of the accretion disc. Without 
attempting to address the issue of including relativistic effects into 
our model in a systematic way, we show here in Figure [8| one of 
our simulated spectra convolved with a relativistic disc model. The 
relativistic blurring routine kdblur (Fabian & Johnstone, priv. 
comm.) makes use of the LAOR kernel describing relativistic ef- 
fects on the spectral shape resulting from emission in an accretion 
disc orbiting a maximally rotating Kerr black hole (Laor 1991). 
Indeed, even for our relatively modest choice of blurring parame- 
ters (emissivity index = 3 and inner radius of the disc at 4.5-Rg), 
most of the discrete emission features in the X-ray energy range 
appear strongly blurred and hardly distinguishable from a curved 
continuum. A more extensive analysis of the strong relativistic cor- 
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Figure 8. The spectrum obtained for ^ = 300 and tb = 3 (red), and the 
same spectrum convolved with a relativistic disc model (in blue, kdblur 
with Ri„ = 4.5-Rg, Rout = 100-Rg, emissivity index 3, inclination of 30 
degrees). The top panel shows the broad band ratio spectra (1 eV-1 MeV). 
The bottom panel focuses on the 0. 1-30 keV range where most X-ray in- 
struments operate. 

rections to our inhomogeneous model is beyond the scope of the 
present paper and will be presented elsewhere. 

In terms of global spectral energy distributions, the general 
trend follows the expectations of the simple analytic two-zone 
model presented in Section|2| Strongly inhomogeneous discs with 
large tb are dominated by reflected/reprocessed UV emission. 
To quantify this, we have calculated the Qox parameter, i.e. the 
slope of spectrum between 2500 A= 5 eV and 2 keV: Oox ~ 
—0.3838 log(-F2kcv/f250o) and the resulting values are shown in 
each panel of figure|6| We see a trend of larger Qox for more com- 
plex, reflection dominated spectra, as indeed observed (Gallo et al. 
2005). At face value, Oox depends quite sensitively on tb, thus on 
the degree of inhomogeneity, and can in principle represent a very 
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Table 1. Spectral parameters obtained from the fit of the simulated spectra 
that illustrate the effects of varying the heating inhomogeneity parameter h 
and the Thomson optical depth of the hot phase . In all the simulations 
tb = 1, 5 = 300. Vp and R are the PEXRAV photon index and reflection 
amplitude, while EW is the equivalent width in eV of the Gaussian line. 
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useful way to tie down the general properties of the observed emis- 
sion. Our model would then predict that the larger Qox, i.e. the X- 
ray weaker a source is, the larger the reflection fraction should be, 
and in general, the more complex the X-ray spectra would appear. 
It has to be kept in mind, however, that the observed emission at 
2500 A may include a contribution from an outer, colder homo- 
geneous accretion disc. If the extent of the inhomogeneous region 
of the accretion disc coincide with the region where radiation pres- 
sure dominates, the relative contribution to the optical emission at 
25OOA should depend on the accretion rate: the larger the accretion 
rate, the larger is the contribution from the inner inhomogeneous 
part we have discussed here. A full model for the dependence of 
Qox on the source luminosity should then be made based on the 
combination of a inner inhomogeneous disc (Appendix 1X1 and an 
outer standard, gas pressure dominated one. 

As expected, the reflection features are weaker when h or tt 
are larger. To confirm this, we have performed two additional simu- 
lations (not shown in Figure and a gain fitted the emerging spec- 
tra with PEXRAV. The results as shown in TableQ If we double the 
size of the inner heated region (i.e. we take h = 0.2), while keep- 
ing tb = 1 constant, the measured reflection fraction, for example, 
drops from 1.73 to slightly less than unity. 



4 DISCUSSION 

The expected variability properties of the inhomogeneous flow 
were discussed in Section ITTI based on the analysis of our sim- 
ple two-zones toy model. Barring any additional (likely) variability 
associated with an intermittent heating rate, we have shown there 
that a modest 10% variation in the optical depth of the cold clouds 
around a mean value of tb ~ 1-2 can produce large amplitude 
variations (up to 40-50%) in the X-ray power-law flux with much 
smaller (less than 20-25%) variations associated with the reflected 
component (see figure |4}. This is a possible explanation for the 
variability pattern observed in many AGN on short (dynamical) 
timescales, as discussed in more detail in Section 1431 Of course, 
as in the geometrically equivalent light bending model of Miniutti 
& Fabian (2004), also in our model there are different regimes of 
variability, essentially depending on the degree off inhomogeneity 
of the heating rate, h, and on the optical depth of the cold clumps 
Tb: when the latter is large, most of the X-ray emission and of the 
reflected radiation are hidden from view, and large amplitude vari- 
ability is induced in both spectral component. On longer (thermal, 
viscous) timescales, the global geometry may change substantially. 
Large variations of tb would of course correspond to substantial 
spectral variability. In figure|9|we have plotted the ratio of the spec- 
tra obtained for ^ — 300 and tb — 1, 5, 8 to the tb ~ 3 spectrum. 
Apart from the dramatic variation of the lowest absorption edge at 
~ 45 eV, and of the overall slope of the Comptonized X-ray emis- 
sion (getting steeper for large tb), the spectra show relatively little 
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Figure 9. Ratio of the spectra obtained for tb = 1 (green), -7-3=5 (blue) 
and Tb = 8 (red) to the spectrum for tb = 3. The top panel shows the 
broad band ratio spectra (1 eV-1 MeV). Solid lines are for ^ = 300, dashed 
lines for § = 3000. The bottom panel focuses on the 0.1-30 keV range 
where most X-ray instruments operate. In this limited energy band, large 
changes in tb can lead to dramatic changes in luminosity with only weak 
changes of the spectral shape. 



variability despite large flux variations, especially in the X-ray en- 
ergy range (0.1-30 keV). In particular, no strong soft excess below 
1 keV is seen in the ratio spectrum. 

4.1 Comparison with earlier theoretical works 

Ballantyne, Turner & Blaes (2004) and Ballantyne, Turner & 
Young (2005) have studied X-ray reflection from inhomogeneous 
discs. The vertical density structure of the disc they assume is de- 
rived from a two-dimensional, time-dependent numerical radiation 
magneto-hydrodynamical (MHD) calculation. However, their nu- 
merical setup was such that the source of X-ray illuminating ra- 
diation was assumed to be external to the disc (a corona above 
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it), and dynamically decoupled from it. Here we have tried to ad- 
dress the different (and computationally much more challenging) 
situation in which there is no external hot corona (and no exter- 
nal energy source), but the heterogeneous, two-phase structure of 
the flow itself gives rise to all the observed high energy radiation 
spectrum. The spectroscopic studies of Ballantyne and coworkers, 
however, suggest that any spectrum produced by a realistic, time 
varying inhomogeneous accretion disc, will pose serious interpre- 
tative challenges to planned future reverberation mapping observa- 
tions. We believe that the exploration of ever more realistic model 
of inhomogeneous discs, possibly including both the density struc- 
ture from 3-D MHD simulations and the inhomogeneous heating 
and the detailed radiative transfer calculation similar to those we 
have presented here will be an extremely useful tool for the next 
generation of X-ray spectroscopes. 

One of the basic assumptions common to many of the differ- 
ent inhomogeneous solutions for accretion discs (see e.g. Krolik 
1998; Begelman 2002; see also the discussion in Appendix 1X1 is 
that the clouds are not confined by magnetic field, but, on the con- 
trary, that magnetic field lines connect the two phases. This implies 
that, for small enough clouds, each of them will be in pressure equi- 
librium with the surrounding hot phase. Then, the density contrast 
(which is a physical quantity easy to derive from numerical MHD 
simulations) will be approximately given by the ratio of the temper- 
atures of the two phases, which could in principle be inferred from 
accurate spectroscopic studies in both soft (0. 1-2 keV) and hard (5- 
30 keV) X-rays. In general, a more realistic density structure that 
takes also into account pressure equilibrium between the phases in 
the presence of a turbulent magnetic field should be the next step, 
although an extremely complicated one, in the implementation of 
our inhomogeneous model. 

The coupled radiative transfer code we have used in this work 
is probably the only one used so far capable of self-consistently cal- 
culating the radiative equilibrium between the two phases present 
in the innermost region of an accretion disc with such a detail. It is 
however very time consuming, and still limited in its wider appli- 
cability. The setup we have discussed here, for example, assumes 
that all clouds are optically thick, and neglects transmission of ra- 
diation across them. If in more realistic inhomogeneous discs the 
cold phase is distributed among clouds spanning a range of densi- 
ties, it is likely that not only reflection off optically thick clouds, but 
also transmission of radiation should play a role in determining the 
observed spectra. In general, in these cases we can expect, for any 
given value of tb, reduced cooling and harder spectra. Moreover, 
the effects of both total and/or partial absorption (depending on the 
filling factor of the cold clumps and on their geometry) will play 
an important role. This interesting issue is far beyond the scope of 
the present paper, which is intentionally aimed at exploring a very 
specific physical configuration. We just note here that a pure "ab- 
sorption model" to explain the soft excess and the general X-ray 
spectral shapes of QSOs and AGN has been recently proposed by 
Gierliiiski & Done (2004b), and a more general and detailed study 
of the role of absorption has been presented in Chevallier et al. 
(2005). 

Finally, we note that more extreme effect of bulk Comptoniza- 
tion which arises in the optically thick plasma in the vicinity of a 
black hole (Blandford & Payne 1981a,b; Payne & Blandford 1981; 
Niedzwiecki & Zdziarski 2006 and references therein), has been 
ignored. However, as in the inhomogeneous flow we envisage the 
great majority of the mass is concentrated in the cold phase (even 
though most of the volume is filled by the hot one), this effect 
should be safely negligible. 



4.2 Comparison with observations: the case of NLSl 

Narrow Line Seyfert I galaxies, identified by the unusual narrow- 
ness of the broad component of their H/3 lines, are believed to be 
powered by SMBH of relatively small masses, with high accretion 
rates, possibly close to the Eddington level (see e.g. Pounds, Done 
& Osborne 1995; Boiler, Brandt & Fink 1996). This fact already 
makes them very good candidate systems for inhomogeneous ac- 
cretion. Furthermore, NLSI have long been known to be charac- 
terized by extreme properties of their X-ray emission: a strong soft 
excess in the ROSAT soft band (0.1-2.4 keV; Boiler et al. 1996); 
unusually steep X-ray spectra in the hard X-ray band (2-10 keV; 
Brandt, Mathur & Elvis 1997; Leighly 1999b, Vaughan et al. 1999); 
very rapid and large variability (Leighly 1999a). 

Recent XMM detailed spectral studies have revealed more un- 
usual spectral properties, most notably in the form of sharp spec- 
tral drops above 7 keV, the most extreme examples of which are 
found in IH 0707-495 (Boiler et al. 2002) and in IRAS 13224-3809 
(Boiler et al. 2003). These features are sharp and time-variable 
(Gallo et al. 2004). Detailed spectral modeling of these data so far 
suggested that either partial covering (Tanaka et al. 2004; Gallo et 
al. 2004) or ionized reflection dominated discs (Fabian et al. 2002), 
with light-bending effects (Fabian et al. 2004; Crummy et al. 2005) 
can explain the observed features in these as well as in other bright 
Seyfert galaxies (e.g. IH 0419-577, Fabian et al. 2005; MCG-02- 
14-009, Porquet, 2005). Within the partial covering model, vari- 
ability is induced by rapid changes in the covering fraction of the 
absorbers. In the light bending model, instead, variability is essen- 
tially produced by a change in the distance between the X-ray emit- 
ting flares and the reprocessor very close to the black hole. 

Our model falls within the family of models trying to explain 
the complex X-ray spectra of bright AGN with ionized reflection, 
though representing an alternative, specific, physically motivated, 
geometrical setup for the multi-phase structure of the inner regions 
of an accretion flow near the Eddington luminosity. In Section |2| 
we have already discussed the qualitative similarities between the 
emerging spectra expected from our inhomogeneous flow and the 
light bending models. Here we would like to notice a few points 
concerning the partial covering model. 

First, a comprehensive analysis of XMM data of a number of 
NLSI with the partial covering model (Tanaka, Boiler & Gallo 
2005; see also Crummy et al. 2005) highlights the fact that the 
estimated temperature of the soft X-rayAJV emission component 
(the soft excess) is very similar for objects that span more than four 
orders of magnitude in luminosity, as already noticed by Walter & 
Fink (1993); Czemy et al. (2003) and Gierlinski and Done (2004b). 
As pointed out by Ross & Fabian (2005), the emission in the 0.2-2 
keV band due to lines and bremsstrahlung in the hot surface layer 
of the cold medium is a signature of ionized reflection, and can 
appear in real data as blackbody emission, provided that the indi- 
vidual lines and edges are smeared by a large enough amount, so 
that the emerging soft excess will appear featureless (Crummy et 
al. 2005; Chevallier et al. 2005). Secondly, the observed rapid and 
strong variability of many NLS 1 puts tight constraints on the max- 
imum distance the partial coverers can possibly be at. In fact, as 
already discussed in detail in Gallo et al. (2004), within the par- 
tial covering model the absorber must be outflowing, have a high 
column density and a low ionization, which imply a kinetic energy 
outflow almost two orders of magnitude larger than the observed 
luminosity. It is more reasonable to assume that the absorber is in 
fact embedded within the X-ray emitter and shares with it the high 



X-ray spectra of inhomogeneous accretion flows 1 1 



rotational velocity of the inner accretion flow, as we have discussed 
here. 

Naturally, reflection dominated inhomogeneous flows from 
the vicinity of a black hole should produce strongly relativistic 
lines in the soft X-ray range. Previous claims of detections of such 
lines in the XMM spectra of two Seyfert galaxies (Mkn 766 and 
MGC-6-30-15, Branduardi-Raymont et al. 2001; Sako et al. 2003), 
have however been questioned (Ballantyne, Ross & Fabian, 2002; 
Turner et al. 2003a) and it is yet unclear to what extent weak rel- 
ativistic lines can be seen in spectra clearly dominated by (dusty) 
warm absorbers (see Lee et al. 2001; Turner et al. 2001). 

In a recent work, Chevallier et al. (2005) have compared ab- 
sorption (and partial covering) and reflection models with the X-ray 
spectra of PG quasars. They reached the conclusion that absorption 
models seem capable of explaining some of the soft X-ray excesses, 
provided the absorbing clouds span a range of optical thicknesses, 
ionization parameters and covering fractions. On the other hand, 
ionized reflection models necessitate that the continuum emission 
is strongly suppressed, as in our inhomogeneous model or in the 
most extreme versions of the light bending model. 

4.3 The scaling with the black hole mass 

The scaling of the physical quantities most relevant for the observ- 
able characteristics of inhomogeneous accretion flows are shown 
explicitly in Appendix A. As the solution we describe is thought to 
be the final outcome of non-linear MHD instabilities in radiation 
pressure dominated discs, it is hard to predict how the parameters 
describing the inhomogeneity, tb, e and h will scale with mass. 
This is in fact the reason why we have chosen them as free param- 
eters of the model. However, once their values are fixed (within the 
range we have explored here), it is apparent that the hot phase tem- 
perature is almost independent on black hole mass (eq. IA13t ,as it 
was already noticed by Krolik (1998). In general, the properties of 
the hot phase turn out to be very insensitive to the central object 
mass. For example, the ratio of gas to radiation pressure in the hot 
phase does not depend on m, ensuring the analytic treatment pre- 
sented in appendix A remains valid throughout the mass spectrum. 
As for the other relevant physical quantities, both cold cloud tem- 
perature and ionization parameter scale with mass as m^^'^*, while 
rici oc m^^''*: stellar mass black holes, as expected, have denser, 
hotter and more highly ionized clumps. 

If one considered also the evaporation/condensation equi- 
librium between the phases (neglected here), as done by Krolik 
(1998), one would infer that the cold phase covering factor should 
rise towards unity as the black hole mass decreases. This corre- 
spond to very large values of tb, for which the observational dif- 
ferences between an inhomogeneous flow and a more standard op- 
tically thick homogeneous disc become less and less pronounced. 

Finally, there is a more general issue regarding the very ex- 
istence of strong inhomogeneities in the accretion flows close to 
the Eddington level. In fact, the critical accretion rate above which 
standard Shakura-Sunyaev discs are viscously and thermally unsta- 
ble depend both on the nature of the viscosity prescription and on 
the mass of the central accreting object (see Merloni and Nayakshin 
2006). The observational fact that accretion discs in luminous black 
hole binaries in the so-called high state appear rather stable close to 
the Eddington rate, (Gierlihski & Done 2004a), but unstable above 
it (Done, Wardzihski, Gierliriski 2004) has been in fact used by 
Merloni & Nayakshin (2006) to argue for a modified viscosity law, 
based on the results of an extended local stability analysis. How- 
ever, if the nature of such a law is unchanged for supermassive 



black holes, the critical accretion rate should scale as m~^^^ . Thus, 
accretion flows in AGN should be unstable at the same (Eddington 
scaled) accretion rates at which high state BHXRB are observed to 
be stable. 



5 CONCLUSIONS 

We have discussed the expected X-ray spectral and variability prop- 
erties of black hole accretion discs at high luminosity, under the 
hypothesis that radiation pressure dominated discs are subject to 
violent clumping instabilities and, as a result, have a highly inho- 
mogeneous structure. Within this picture, most of the mass is con- 
centrated in (Thomson) optically thick clouds that fill only a small 
fraction of the volume, which is instead occupied by hot, optically 
thin Comptonizing plasma. Such a geometrical configuration was 
first proposed by Fabian et al. (2002) as a mean to explain the com- 
plex spectral properties of the extreme NLSl IH 0707-495. Since 
then, many more bright AGN (both Seyfert and QSO) have been 
modeled with reflection dominated spectra (see e.g. Crummy et al. 
2005). 

Here we have presented a physically motivated model for ac- 
cretion discs near the Eddington luminosity (Appendix |^, and 
studied their radiative output by exploring the space of parame- 
ters that describe the proposed geometry (mean optical depth of 
the cold phase and degree of concentration of the heating source) 
in a simple two-zone analytical model. Moreover, we have used a 
sophisticated radiative transfer code to simulate the radiative equi- 
librium between the two phases in a self-consistent way, so that the 
mutual radiative feedback is fully taken into account. 

The main results of our work are the following: 

• for very inhomogeneous configurations, i.e. those which have 
both a large covering fraction in clouds and very concentrated dis- 
sipation, the emerging radiation can be dominated by ionized re- 
flection; 

• the observed spectra are very complex, with prominent signa- 
tures of ionized reflection dominated by a plethora of UV and soft 
X-ray emission lines; if the blurring due to relativistic motions in 
the immediate vicinity of a black hole is accounted for, the soft X- 
ray emission will appear as a featureless continuum (a soft excess); 

• for a large range of parameters, even when contributing a large 
fraction of the hard X-ray spectrum, the reflection component is 
less variable than the power-law like emission originating from the 
hot Comptonising phase, in agreement with what observed in many 
Narrow Line Seyfert 1 galaxies and bright Seyfert I; 

• the reflection fraction becomes strongly variable only at very 
low flux levels of the hard X-ray continuum, corresponding to a 
configuration in which both the primary Comptonized continuum 
and the high energy reflection emission are essentially hidden from 
the observer's view. 

The fact that many of these properties are common to the 
"light bending" model (Miniutti & Fabian 2004) is not surpris- 
ing. In fact, for a two-phase system in radiative equilibrium, the 
main spectral properties are determined by the radiative feedback 
between the hot and the cold phases. Such feedback is strongly de- 
pendent on the geometry (and on the topology) of the two phases 
and, in particular, on the sky covering fraction of the cold phase 
as seen by the hot, Comptonising medium. Reflection dominated 
spectra are expected when the cold phase intercepts most of the 
photons coming from the hot phase. 

In more general terms, one can define "open" any geometry in 
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which the hot, X-ray emitting plasma is photon starved. Such a ge- 
ometry will produce hard X-ray spectra, little soft thermal emission 
and weak reflection component. On the other hand, a "closed" ge- 
ometry will instead correspond to a very large covering fraction of 
the cold phase, with associated strong soft emission, softer spectra 
and strong reflection fraction. Such a geometry corresponds, obser- 
vationally, to the highly accreting black holes and is obtained via 
general relativistic (GR) effects in the light bending model, while 
is a result of the clumpy and inhomogeneous nature of the inner 
unstable, radiation pressure dominated part of the disc in the model 
we have presented here. 

In the end, it seems likely that a complete model for the in- 
nermost region of luminous accretion discs around black holes will 
involve a two-phase inhomogeneous and turbulent flow moving at 
relativistic speed, so that the emerging spectrum will be determined 
by ionized reflection off the cold phase, inverse Compton scattering 
in the hot phase and possibly warm absorption from an accompany- 
ing outflow, all distorted by general relativistic effects. The model 
we have explored here, thus, represents a step towards our under- 
standing of such a challenging complexity. 
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APPENDIX A: A FULL INHOMOGENEOUS ACCRETION 
DISC MODEL 

We derive here a full accretion disc model in which the inhomoge- 
neous structure is introduced through the parameters tb, tt, 6, e, 
h, as discussed in the paper. We will follow closely the approach 
of Krolik (1998) that studied analytically the overall equilibrium 
of a clumpy two-phase accretion disc under the hypothesis that the 
cold clumps are magnetically connected to the hot phase. This con- 
figuration turns out to be both thermally and viscously stable. Its 
energy flow is more complicated than that of standard uniform disc 
as exemplified by the flow chart of Fig.l of Krolik (1999). The 
hot phase heat comes either from magnetic field reconnection and 
thermalization of cloud random motions (either due to drag or in- 
elastic collisions), while the cold phase is heated either through 
cloud-cloud collisions and through reprocessing of the hot phase 
radiation. In order to allow a simplified treatment of the radiative 
coupling between the phases, both in our toy model (Section|2j and 
in the radiative transfer simulations (Section|3j we have neglected 
internal dissipation in the clouds due to their inelastic collisions, 
and focussed only on their heating by absorption of radiation emit- 
ted in the hot phase. We will keep this assumption here, equivalent 
to setting to zero the term ac in Krolik (1998). 

We begin by expressing the angular momentum conservation 
equation in terms of the cold clumps velocity dispersion a^i- Here 
we adopt the customary a viscosity prescription of Shakura & 
Sunyaev (1973), independently of the details on how exactly the 
stress is produced (see Krolik 1999 for a discussion). In the fol- 
lowing, we will indicate the accretion disc radius in units of the 
Schwarzschild radius as r, and the accretion rate in terms of the 
dimensionless quantity m — eradMc^/I/Bdd with the radiative ef- 
ficiency erad = 1/12 fixed to its Newtonian value, as appropriate 
for the inner boundary condition we have adopted. We have: 

«<7ci = — (Al) 

where S7 is the Keplerian angular velocity, E — 2pH and the func- 
tion J(r) = (1 — fin/r), with rin = 3, describes the Newtonian 
approximation of the no-torque at the inner boundary condition for 
a disc around a Schwarzschild black hole. The velocity dispersion 
of the clouds can be expressed in terms of the hot phase sound 
speed, Cs.h'- 



3.2 X 10^e]/^M, 



(A2) 



where is the temperature of the hot phase in units of m^c^ /k. As 
discussed in Krolik (1998), if the mean separation between the 
clumps is of the order of the disc scaleheight, the difference in the 
potential between neighboring clumps is then big enough to stir 
them aX. M ~ 1. In fact, within our simple geometrical model of 
Section|2|and for small enough clumps, we can estimate the mean 
clump separation in units of iif as ~ 0.2Tg^^^ . 



The angular momentum conservation equation <A1> can then 
be cast as a constraint on the total (cold plus hot phases) average 
Thomson optical depth of the inhomogeneous disc, rtot = kE/2, 
with K electron scattering opacity, i.e. the optical depth the disc 
would have if it was homogeneous: 

rtot = 2.6 X 10^Q~V~='''^mJ(r)7W"^eil . (A3) 

Let us now discuss the hydrostatic equilibrium of the disc. For 
the hot phase, assuming that gas pressure dominates over radiation 
pressure, the disc scaleheight is given by: 

Ji- ~ 1.4 X lO^Vvr^/^ei/i^ . 



(A4) 



The cold clumps, obviously, will be distributed over a scaleheight 
of the order MH. 

We can now verify that indeed our assumption about the pres- 
sure dominance in the hot phase is justified. The gas pressure is 
given by: 

Pgas.h = 2phfcrhot 1.9 X lO^mf V-^/'ei^iVT dyne, (A5) 
while for the radiation pressure we have 

-Prad.h = rTQ/(3c) ~ 1.8x 10** m7V"^m_iJ(r)rT dyne,(A6) 

where Q is the disk dissipation rate per unit area. Their ratio is then 
given by: 

-frad,h 



5.8r 



-3/2^ 



iJ{r)QZ^^ < O.llm-ieii''^ , (A7) 



1/2 



where the inequality in the above expression comes from taking the 
maximum of the radial function r~^/V(r) (at r = 16/3 ~ 5.33). 
So our assumption is justified for sub-Eddington disks. 

The radiative properties of the system will be in general a 
function of the density, temperature, total number and geometry 
of the cold clouds. We can estimate the cold clump temperature un- 
der the zeroth order assumption that they will re-radiate thermally 
all the radiation coming from the hot phase. Consequently, we can 
write Q/{ch^) — AasThh/c, where ob is the Stefan-Boltzmann 
constant, and we have allowed for inhomogeneous heating by in- 
troducing the parameter h, the size of the heated region in units of 
the disc scaleheight (see Section|2j- Thus we get, for the cold cloud 
temperature: 

fcTbb ~ 1.1 X 10^m-'''V-^''''[(l-a)/0.64]'/* 

X [rn-iJ{r)]^'^hZ\'^ eSf . (A8) 

Then, assuming pressure balance between the clouds and the hot 
ambient phase we get, for the cloud density: 



"-cl 



1.0 X 10'^m-^''V-^/''[(l - a)m--,J{r)]-^'* 



, 1/2 ^1/2 -3 

X h_-^ ttQ_i cm . 



(A9) 



Neglecting the weak dependence on the albedo a, the typical ion- 
ization parameter, in units of erg cm s~^, is given by 

C~ 1.3xlO"m-^/V-'*/*[m_iJ(r)]^''^ft:'^/V-^e:y' .(AlO) 

For the fiducial values of the parameters we are discussing here, ^ 
attains a maximum value of about 10"^, close to the radius of max- 
imal dissipation in the disc. Also as expected, the strongest depen- 
dence is on the geometrical factor h, as the ionization state of the 
cloud surfaces will be most sensitive to the heating inhomogeneity. 

The cold clouds filling factor is given by / = feTs, and the 
total density in the disc is then 

Tbb 4 

Wtot = fnci + {l-f)nh = nci[f+{l-f)— — ] ~ -ersrid .(AH) 

J hot 3 
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It is now possible to close the disc structure equations and find a 
relationship between the model free parameters (rs, tt, 6, e, h) 
and the physical parameters of the standard accretion disc solutions 
(mass, accretion rate, viscosity, etc.). In order to do so, we equate 
the total disc optical depth, dominated by the cold phase, arritotH 
to the expression in equation <A3> . We obtain, finally: 

TBTTe'iie^ih]/^ ^ 22r-^^''[m.,J{r)f^^m-'^^a-^M-\Al2) 

So far, we have neglected the condensation/evaporation equi- 
librium between the cold and the hot phase. If we introduce the 
relationship between the typical cloud size and the hot phase tem- 
perature and density, as derived by Krolik (1998, eq.5), we obtain 
the following expression for the hot phase temperature: 

e ~ 0.36Q-^''"m-^/"r-'-*/"[m_iJ(r)]^/" 

X M'^^'r-^/'h]/^ , (A13) 

almost independent on black hole mass and only weakly dependent 
on accretion rate. 

To conclude, it is interesting to notice that we have obtained 
a self-consistent solution for the inner radii (i.e. those where radia- 
tion pressure is supposed to dominate) of supermassive black holes 
accreting at a few tenths of the Eddington rate for our fiducial pa- 
rameters, i.e. both Tt and tb of the order unity and h, e, B all of 
the order 10^^ . In both Sections|5|and0we have explored in detail 
the radiative equilibrium and the expected spectral and variability 
signature of an accretion disc in which all the above parameters 
are either fixed at, or allowed to vary around, these fiducial num- 
bers. The fact that they indeed correspond to the right values for a 
full accretion disc model as we have outlined here is indeed a very 
encouraging result. 



APPENDIX B: ANALYTICAL ESTIMATES OF THE 
EMERGING LUMINOSITIES FOR THE TWO-ZONE 
MODEL 

For each zone, the escaping luminosity is related to the internal en- 
ergy through the escape probability (Lightman & Zdziarski 1987): 

P={L/U)(H/c), (Bl) 

which depends on the geometry, sources distribution, optical depth 
and energy. We consistently indicate as Pc"', Pvy'' ,1^1^ , the es- 
cape probability for Comptonized, soft and reflected radiation, with 
j = i, o for the inner and outer zone, respectively, and will use the 
analytical approximation to P given in the Appendix of MC02. For 
the inner zone we assume for the escape probability the analytic 
formula appropriate for a spherical geometry with central isotropic 
injection (eq. A3 of MC02). For the outer zone, a central isotropic 
injection in spherical geometry is appropriate as long as ft <C 1, 
while for h — > 1 a slab geometry should be more appropriate. 
Thus, we employ an interpolating expression for the intermediate 
cases. 

Also, we should include a term that depends on the average 
fraction K of the emerging luminosity from the outer zone that 
reenters the inner one. This is computed as follows: consider a shell 
of width dr located at a distance r from the center, with r > h. The 
solid angle (divided by 2n) subtended by the inner zone as seen 
from the shell is K(r) = Afi/27r = 1 - y/l - (h/r)'^. Then we 
can define 

K = K{r)dr = 1 - J/(l - h), (B2) 



where the integral is evaluated to 

Vl-/t2+/iarcsin/i-ft(7r/2). (B3) 

As we assume uniform density for the hot and cold phase, we 
have that the optical depths in the two zones are simply given by: 
= h-rr; = (1 — h)TT, and analogously for re. 

Inner Zone: The soft luminosity in the hot phase is produced by 
the clouds through absorption and reprocessing of the Comptonized 
and reflected radiation and disappears through escape or Comp- 
tonization, additionally, a fraction of the soft luminosity emitted by 
the outer zone reenters the inner zone. Thus, the radiative equilib- 
rium balance for U\jy reads: 

(Puv' + rT)(7uv = T^{1 - a)Uh + r^{l - aR)UR 

+ (if)Puv°f/5v, (B4) 

where a and a_R are the clouds energy and angle integrated albedo 
for a Comptonized and a reflection spectrum, respectively. 

Similarly, the reflected radiation is formed through reflection of 
the Comptonized and reflected radiation on the clouds, and disap- 
pears via escape, Comptonization in the hot plasma and absorption 
by the clouds: 

[Pn +r^ + riil - aB)]Un = r^aUh + KP^Ur. (B5) 

The fraction of the internal energy in the form of Comptonized 
radiation has as main source the power dissipated in the hot plasma 
but also the Comptonized soft and reflected radiation; the sink term 
is due to escape, absorption and reflection by the clouds: 

{Pc + rl,)Uh = I^hi^/c + T^Uh^ + rirU^B + KPc°U^. (B6) 

Outer Zone: Similarly, we can write the equations for the radia- 
tive equilibrium in the outer zone, considering that all the luminosi- 
ties emerging from the inner zone have to enter the outer zone: 

{PvN° +r^)U^Y = T^il - a)Uc + r^O^ - aR)UR 

+ PuvUhv, (B7) 

[Pr +r^ + rg(l - aR)]U°R = rgat/g + PrU'r, (B8) 

(Pc° + T^Wc = r^U?,v + T^U°R + PcUh- (B9) 

The above system of six equation in the six unknown ((7/, 
with j = i, o and i — C, R, UV) can be solved analytically, and 
the three escaping luminosities finally computed through eq. <BH . 

An estimate for the observed reflection coefficient can be di- 
rectly obtained from: 

(BIO) 

where Lx is the Comptonized luminosity above 1 keV. For the net 
albedo (accounting for multiple reflections), Onct, we use the ex- 
pression in eq.j4}, derived from a calibration with the numerical 
simulations, as discussed in Section|2| This formula corresponds to 
an angle averaged reference slab-reflection spectrum. Actually, the 
reflection coefficient also depends on the assumed inclination angle 
for the infinite slab model (see e.g. PEXRAV model in XSPEC, 
Magdziarz & Zdziarski 1995). Equation <B10t could be corrected 
for this, for example by dividing it by the angular factor given by 
equation 2 of Ghisellini, Haardt & Matt (1994). However, for the 
relatively low inclination angles usually assumed in spectral fits, 
the correction is small and, for simplicity, we will neglect it. 
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